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We investigate the propagation of a wave-packet in the c/> 4 model. We solve the 
time-dependent equation of motion for two distinct initial conditions: The wave- 
packet in a trivial vacuum background and in the background of the kink soliton solu- 
tion. We extract the scattering matrix from the wave-packet in the kink background 
at very late times and compare it with the result from static potential scattering in 
the small amplitude approximation. We vary the size of the initial wave-packet to 
CN . identify non-linear effects as, for example, the replacement of the center of the kink. 
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INTRODUCTION 



Oh ' Non-linear models with soliton solutions [1] possess a wide range of applications in 

physics. To present an incomplete list, this range covers cosmology 0, IH, particle and 
nuclear physics @, HJ , as well as condensed matter physics j6[ . The kink-model in one time 
and one space dimension for the real scalar field 4> 

^ . is a prototype model with a topological soliton solution based on spontaneous symmetry 
^ | breaking. Any feature observed in this model is likely to be relevant for the above mentioned 
applications. 

In the Lagrangian, eq. (pQ) M is the mass parameter and A is the coupling constant. The 
soliton solution is the so-called kink 

X ; 4>k{x) = -^Ltanh ( — x J . (2) 



: V2A V 2 

It depends on the spatial coordinate but is time-independent and interpolates between 
the degenerate vacuum configurations 0o = ^^fp that arise from spontaneous symmetry 
breaking. Any deviation from either of these two values at spatial infinity causes the total 
energy 

dxe(x,t) with e(x,t) = - (0 2 (x,t) + /2 (x,t)J + - f 2 (x,t) - — 

(3) 

to diverge. Here we have introduced the partial derivatives as <fi = dcj)/dt and 0' = dcj)/dx. 
For later use, we have made the time dependent energy density, e(x,t) explicit. Since such 
divergences may not occur when continuously converting one configuration into another, 
any smooth deformation cannot alter the field configuration at spatial infinity. In particular 
it conserves the difference </>(— oo,t) — <f>(oo,t). This difference may thus be viewed as a 
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topological charge. With proper normalization, the kink carries topological charge one, 
while the vacuum configuration t) = 0o has charge zero. 

Commonly small amplitude fluctuations about the kink, <f)(x,t) = 4>k( x ) + v( x ^) are 
introduced to determine scattering data. The small amplitude wave-function r](x,t) obeys 
a Schrodinger type equation with a space-dependent potential induced by the kink. In 
this small amplitude approximation quadratic and higher orders in r)(x, t) are omitted when 
substituted into C At spatial infinity, i.e. far away from the kink, this Schrodinger equation 
is a Klein-Gordon equation with mass M. Application of standard methods [l| to this 
Schrodinger problem yield the phase shift 1 

„ n „ / 3Mk \ 

6(h) = 2arctan ) • (4) 

The generalization to three spatial dimensions is that of meson-baryon scattering [3]. The 
kink 0k represents the baryon and the fluctuations 77 are mesons 

In this paper we show how the same result for the phase shift can be obtained from 
the physically motivated set-up of a wave-packet interacting with the kink. This exactly 
matches the scenario that a meson scatters off a baryon. In particular, we construct the 
wave-packet such that at spatial infinity the field equals either of the possible vacuum con- 
figurations for any finite time t. In this set-up we will numerically solve the time dependent 
equation of motion 

00, t) = <p"(x, t) — X U 2 (x, t) - - </>(x, t) (5) 

with appropriate initial and boundary conditions for (f>(x,t). The total energy, eq. (j3J) of 
configurations that obey this equation of motion stays constant in time. 

Numerical simulations of the time-dependent solutions in the kink model have been 
around since quite a while. For example, the transition time of a wave-packet through the 
kink was estimated in ref. This occurred in the greater context of modeling phase tran- 
sitions. The resonance structure of the kink-antikink interaction was numerically studied 



in ref. [lOj. Lattice studies of the kink model were reported in ref. [Ill, [12j as a simula- 
tion of molecular dynamics. These numerical studies focus on kink-antikink configurations 
which belong to the topological trivial sector. More recently this sector has un derg one 



numerical investigation in the context of pair-creation [13[, electro-weak oscillons [14|, [15 



and bounces [16|, |17] • Numerical simulations within the 2 + 1 dimensional generalization at 



non-zero temperature allowed to study the effect of fast quenches to model resonant nucle- 



ation [18[. This scenario may have applications in condensed matter and cosmology [19 
Within super-symmetric extensions of the kink model, numerical simulations have even been 
used to study brane world collisions [20-22||. In the unit topological charge sector numeri- 



cal simulations have been employed to investigate the effect of impurities in the kink [23 
However, no study in the sector of unit topological charge is known to us that attempts 
to extract the scattering phase shift from the time-dependent differential equation or to 
identify non-linear effects in this sector by going beyond the small amplitude approximation 
for the scattering process. 



Since the induced potential is rencctionless, the two eigenphase shifts are identical. 
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This short paper is organized as follows. In the next section we describe the set-up 
for the wave-packet, i.e. we specify the initial and boundary conditions to solve the time 
dependent equation of motion In principle, a single integration of this equation can, 
after Fourier transformation, provide information on scattering data for all momenta. In 
section III we will discuss the numerical simulation of the time-dependent configuration and 
describe how the scattering phase shift can be extracted. We will observe structural changes 
of the interaction pattern as the amplitude of the initial wave-packet increases. We conclude 
in section IV. 



II. WAVE PACKET 

We will consider wave-packets in the two sectors of vanishing and of unit topological 
charges. Their structures are generic and we merely have to specify the background fields 
for either of the two cases. 

As initial condition for the wave-packet we consider a linear combination of plane waves 
that satisfy the Klein-Gordon dispersion relation, u k = \/k 2 + M 2 . Hence at t = we 
parameterize 

/oo 
dkA(k)e ikx . (6) 
-oo 

At t — the dispersion relation enters only via the velocity of the initial wave-packet 

/oo 
dku k A(k)e ikx . (7) 
-oo 

In the context of our numerical simulations we assume the spectral function in momentum 
space to be of Gaufiian shape 

A(k) =a e ~* , (8) 

where ao is the amplitude, ko is the average momentum, and a k is the width of the distri- 
bution. Though the momentum integral in eq. ([6]) can be straightforwardly computed, we 
refrain from displaying it here because its analog in eq. (J7J) has no closed representation. We 
compute either of the integrals numerically. We remark that r] wp (x) may exhibit oscillations 
since it is a complex Gaufiian distribution. This is especially the case when a k is not partic- 
ularly large. The distribution is centered around x = and fj wp (x) increases with fc , which 
measures the average velocity with which the wave-packet propagates. 

As the equation of motion is second order in time, we have completely determined the 
initial conditions for the wave-packet in eqs. ([6]) and ((7j). We next super-impose it with 
various background configurations that solve the static equation of motion. Such combina- 
tions serve as initial configurations to integrate the equation of motion (jSj). By construction, 
the initial wave-packet is localized and since the total energy is conserved, the topological 
charge of this super-imposed configuration will be that of the background. 



i 



A. Pure wave packet 



For the pure wave-packet we merely consider the above defined distribution around the 
trivial vacuum at early times 



0(x,O) 



M 



+ Tj 



wp 



X 



and 



<90(x, t) 



Of 



t=Q 



(9) 



and feed it into the equation of motion ([5]). Since the velocity fi wp (x) is characterized by the 
dispersion relation extracted from eq. ([5]), the superposition 



(j)(x,t) 



M 



+ 



dkA(k)e 



i(kx—ui^t) 



(10) 



is an approximate solution as long as ao is small enough to neglect 0(7]^) terms in the 
equation of motion. This omission defines the small amplitude approximation. As we 
increase ao non-linear effects emerge and the solution to the differential equation will no 
longer be a superposition of plane waves. 

We may converse this line of argument to study non-linear effects. Assume <f)(x,t) to 
be the (numerical) solution to the equation of motion that emerges from the above defined 
initial condition. Then the deviation of the Fourier transform 



4>t s {k) 



Q^ktf 



dx 
2^ 



-ikx 



<j>(x,tf) 



M 



(11) 



from unity measures non-linear effects for times tf ^> 0. Obviously, a single integration of 
the equation of motion in coordinate space provides information about the full momentum 
space. 



B. Wave packet in kink background 

For the wave-packet with unit topological charge in the kink background we consider the 
initial configuration 

(f)(x, 0) = <p K (x - x ) + rj wp (x) (12) 

where Xq is the position of the center of the kink soliton solution. It must be taken large 
enough to avoid any overlap between the kink and the wave-packet at t = if we want to 
discuss the scattering problem. In order for scattering to occur, the signs of Xq and ko must 
coincide. Otherwise the wave-packet will propagate away from the kink 2 . Since the kink is 
static, the initial velocity is as in eq. (J7J). Any non-zero velocity of the kink can eventually 
be compensated by an appropriate Lorentz transformation [lj], that modifies the details but 
not the structure of the wave-packet. 

Again, we can give an analytical expression for the solution to this initial condition, 
provided we may omit 0(r]^ p ) terms according to the small amplitude approximation 

<f>(x,t) =<f> K (x-x )+riW{x,t), (13) 



2 The wave-packet contains components with negative momenta. They do not participate in scattering for 
our choice xq > 0. 
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where 

/oo 
dkA(k)exp[i(kx-u k t + S(k))] . (14) 
-oo 

As for the pure wave-packet, our analysis will be converse. We prescribe the spectral 
function A(k) as in eq. (ED together with the dispersion relation associated with the Klein- 
Gordon equation. This enters the initial configurations 4>(x, 0) and <p(x, 0) and we utilize 
the equation of motion, eq. (jSJ) to find the time-dependent configuration, <p(x,t). We then 
consider very late times tf 3> at which the wave-packet has completely penetrated the 
kink and the two structures are again well separated and can be individually identified. This 
defines the late-time wave-packet type solution 

VwJ(x, tf) = 4>{x, t f ) - <p K (x - x ) . (15) 

Its (inverse) Fourier transform should be compared with the small-amplitude solution, 
eq. flU 

/oo J 
Q 2^ (X ' tf) = m 6XP [l {6{k) - Uktf)] + °^ ■ (16) 

That is, from the numerical solution to the equation of motion we should be able to extract 
the phase shift 

J^tf rco dx 



C 



iS(k) _ f / ^± -ikx (5) 



As long as t](x,t) satisfies the criteria for a small amplitude fluctuation, the dependence on 
tf cancels on the right hand side. A main purpose of the present investigation is to compare 
the numerical result, eq. ( 117]) . for the components that participate in scattering, with the 
result from small amplitude approximation in eq. fl4]). 



III. NUMERICAL RESULTS 

Appropriate scaling of the coordinates and the field 

^ (18, 

allows us to completely absorb the model parameters. Hence their actual values are of minor 
relevance and all results are genuine. In this section we will quote all numerical results in 
terms of the dimensionless quantities on the right hand side of eq. ffl8|) . In these units the 
vacuum solutions are at </> = ±1 and the small amplitude fluctuations have mass a/2. 

The numerical treatment starts by defining an equi-distant grid with spacing h in coor- 
dinate space. This establishes an interval on the x-axis that we take to be finite but much 
larger than the extension of the wave-packet and the kink. Then we employ a fourth order 
Runge-Kutta algorithm together with an adaptive step size control to solve the equation 
of motion (j3J). The latter considerably slows down the numerical computation as the am- 
plitude a is increased. Earlier attempts using a simple Euler algorithm failed to produce 
acceptable accuracy. The equation of motion propagates the configuration in time. At each 
time step (as well as at the auxiliary intermediate steps required by the Runge-Kutta algo- 
rithm) we compute the (second) spatial derivative of the configuration that occurs on the 
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FIG. 1: (Color online) Time snapshots of the energy density e(x, t) (in M 4 /8A 2 ) of the wave-packet 
for real initial conditions. We have used ko = 4 and <Jk = 2. Left panel: a = 0.05, right panel: 
a = 0.15. The numbers next to the lines refer to the time variable. Note the different scales on the 
ordinate. 



right-hand-side of the equation of motion with an 0(h A ) accuracy. To this end, the config- 
uration is assumed to vanish at points outside the considered interval in coordinate space. 
This corresponds to the boundary condition that no flux penetrates outside this interval. 
As a consequence thereof, the wave-packet bounces at the spatial boundaries after very long 
times. This is, of course, not physical but merely a finite size effect and we have to terminate 
the simulation at late times when this phenomena becomes visible. 

A major criterion to accept the numerical solution is that the total energy, eq. ([3]) stays 
constant in time (at the order of the desired numerical accuracy). The complex wave-packet 
initial condition implies the total energy and the energy density to be complex as well. The 
investigation of the physical energy density, e(x,t) in eq. ([3]) hence requires to also solve the 
equation of motion with the real initial condition 

/oo /*oo 
dk A{k) cos(fcx) and Vr( x ) = / dkuo k A(k) sin(fcx) (19) 
-oo J — oo 

for the wave-packet. In figures ([T|) and (T2]) we display the time evolution of the subtracted 
energy density 

e(x,t) = e(x,t) -e bg (x) (20) 

for this initial condition. To single out the wave-packet contribution, we have subtracted the 
energy density associated with the static background. For the pure wave-packet this is zero 

but with the kink background we have £bg(%) — \ 1 — tanh 2 ^^^Jp^ in the dimensionless 

units of eq. ()18p . In the context of the small amplitude approximation the kink is assumed 
to be infinitely massive and thus does not change its location during the interaction 3 . We 
observe an interesting effect at the center of the kink. Once the wave-packet has passed by, 
a residual deformation of the energy density remains. This effect persists even for very long 



3 In the three dimensional scenario of meson baryon scattering this resembles the large Nq picture in which 
the baryon is O(Nc) heavier than the meson. 
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FIG. 2: (Color online) Same as figure Q] with the kink background. 



times and becomes stable. Also, it is the more pronounced the larger the amplitude of the 
wave-packet is. We will discuss a potential explanation for this effect in subsection B. 

The energy density associated with the wave-packet in particular suggests to discuss its 
spread by first defining the normalized expectation values 



(x»> 



/ dx x n e(x,t) 
J dxe(x,t) 



This enables the computation of the (squared) standard deviation 



a' 



(x 2 ) - ( x y 



(21) 



(22) 



as a direct measure for the width of the wave-packet. The time-dependence of the position 
of the center (x) is essentially unaffected by the kink as the data in table UJshow. Furthermore 



its velocity agrees with what is expected for the wave-packet as 



fco 



0.94 for k = 4. 



We find the total energy stored in the wave-packet to be 0.52 and 4.65 for ao = 0.05 and 
ao = 0.15, respectively. The latter corresponds to more than three time the mass of a Klein- 
Gordon particle. Non-linear effects for (x), i.e. its dependence on ao, are only marginal. We 
compare the results for a from different values of the amplitude ao i n both cases, with and 
without the kink background, in table [Til Obviously the kink background causes a significant 
increase of the spread of the wave-packet. Closer inspection shows that this manifests itself 
mainly after the interaction between the wave-packet and the kink, that is for t > 10 while 
even up to the time of the interaction no significant difference between the cases with and 
without the kink is observed. We assert the strong increase of a after the interaction to the 
emergence of the above mentioned structure around x ~ Xq rather than to a direct spread of 
the wave-packet. Indeed, the comparison between figures [T] and [2] does not indicate a severe 
increase of the spread. 



A. Propagation of pure wave packet 

We now return to the complex valued initial wave-packet. We first consider the pure 
wave-packet, eq. (Q. The simulation of eq. f TTTj) will provide information about the numerical 
accuracy that we can expect when attempting to extract the phase shift at a later stage. 
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10 


-0.56 


-0.58 
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10 


0.77 


0.95 


0.80 0.97 


15 


4.15 


4.13 


4.17 4.17 


15 


0.84 


1.16 


1.65 1.81 


20 


8.87 


8.83 


8.89 8.90 


20 


0.91 


1.38 


2.25 2.50 



TABLE I: The central position (x) of 
the wave-packet as a function of time 
for k$ = A and at = 2. 



TABLE II: The normalized standard 
deviation, a, as a function of time. Pa- 
rameters are as for table HI 





FIG. 3: (Color online) Time snapshots of the real part of the configuration 4>(x,t) — M/y2A for 
the initial condition, eq. (|12p . Furthermore we have used k$ = 4 and = 2 to characterize the 
wave-packet. Left panel oq = 0.05, right panel ao = 0.15. Note the different scales on the ordinate. 



Figure |3] shows the real parts of the numerical solution to the equation of motion for 
two different values of the amplitude a Q . The imaginary part behaves similarly, just phase 
shifted by tt/2. Note also that in addition to the expected spread, the number of (visible) 
oscillations contained within the wave-packet increases with time. 

From the comparison of the two graphs in figure [3] no significant dependence on the 
initial amplitude of the form of the wave-packet can be deduced. This confirms the results 
for the standard deviation listed in table [TTJ This absence of significant non-linear effects 
is somewhat surprising. Rather they were expected as a signal for particle production, as 
sufficient energy is available. 

In figure H] we show the inverse Fourier transform denned in eq. (fTTj) . As indicated above 
its deviation from unity provides insight in the numerical accuracy that we can achieve, 
at least for small ao- At low momenta deficiencies arise because these contributions have 
left the bulk of the wave-packet. Furthermore small errors at low and large momenta are 
amplified as 1/A(k) is large in these regimes. The deficiencies at small momenta are obviously 
independent of ao- This is not the case for those at large momenta, even though their 
structures are very similar. In this momentum regime they originate from the inverse Fourier 
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k 



FIG. 4: The inverse Fourier transform of eq. (jlip for ao = 0.05 (left) and ao = 0.20 (right). The 
dashed vertical lines separate regimes with different numerical parameters to improve the numerical 
accuracy. 

transform requiring a very fine grid, i.e. small h, in coordinate space. Unfortunately this 
comes with a heavy computational cost and it occurs appropriate to divide the momentum 
axis into subintervals that are treated with different numerical parameters. This not only 
concerns the parameters for numerically solving the equation of motion (e.g. the step size h) 
but also the detailed structure of the wave-packet that we characterize by k$ and o^. At the 
interface of these subintervals the results extracted from the different numerical treatments 
match. This is also indicated in figure HI At this stage it is difficult to judge whether the 
deficiencies at large k are of numerical origin or signals of the non-linear dynamics. 

B. Propagation in kink background 

In figure we present solutions to the numerical integration of the equation of motion 
at different times for small and moderate initial amplitudes of the wave-packet. At all 
times, the deviation from the kink is confined within a spatial regime. We clearly identify 
the interaction process when the wave-packet climbs up the kink. A wave-packet with 
a small amplitude essentially retains its shape after the interaction, the only effect being 
characterized by the phase shift that we will extract later. Surprisingly, even for moderate 
amplitudes the shape of the wave-packet does not change significantly with time. 

Figure \5\ has the potential to explain why the energy density e developed some structure 
around x = Xq after the interaction with the kink: The kink acquires a displacement d. 
This displacement increases with the initial amplitude a . We quantify this increase by 
identifying the point d at which </>(d, t) — in the vicinity of Xq for very late times. We 
list d for various values of ao < 0.2 in table IIHI In agreement and support for the small 
amplitude approximation, the displacement vanishes with the amplitude. For moderate ao it 
saturates after a while and stays constant; at least at the order of our numerical accuracy 4 . 
As ao increases further, the displacement changes sign and slowly grows with time. This is 
the onset of a novel non-linear effect that we will discuss further in subsection IIIID1 Our 
numerical simulations suggest that d does not depend on k . 



4 On the overall scale this displacement is an effect at the order of a fraction of a percent. 
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FIG. 5: (Color online) Time snapshots of the real part of the configuration </>(x,t) for the initial 
condition, eq. (fl~3l) . Left panel ao = 0.05, right panel ao = 0.15. 





t=100 


t=150 


t=200 


a 


0.05 0.10 0.15 0.20 


0.05 0.10 0.15 0.20 


0.05 0.10 0.15 0.20 


d 


-0.07 -0.40 -0.80 1.32 


-0.07 -0.39 -1.13 2.80 


-0.07 -0.40 -1.09 3.05 



TABLE III: Displacement, d measured relative to xq, of the kink as a function of the initial 
amplitude of the wave-packet, ao- 



For negative d the numerically observed structure (see figure [2]) of e(x, t) in the vicinity 
of Xq is reasonably reproduced by e^x — d) — 6b g (x), i.e the corresponding shift of the 
background energy density. The larger the time interval the better is the agreement with this 
analytic expression. To separate the displacement effect and to focus on the time-evolution 
of the wave-packet we repeat the calculation of eqs. (12Tj) and with the lower boundaries 
of the integrals taken between the kink and wave-packet as we can separate these structures 
unambiguously. The results are shown in table HVl The comparison with the time-evolution 
of the wave-packet in the trivial background reveals that the interaction with the kink has 
only a minor influence on the spread of the wave-packet once the displacement of the kink is 
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178.7 178.4 


5.09 9.92 


179.3 180.2 


4.94 11.40 



TABLE IV: Center and spread of the wave-packet in the kink background for very late times. We 
compare the two cases with the vacuum and the kink backgrounds. In either case we have used 
ko = 4 and = 2 with real initial conditions. 
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properly accounted for. In both cases the spread increases with the amplitude of the wave- 
packet. This summarizes the main features of the wave-packet for small and moderate 
amplitudes. For a > 0.2 we observe a different behavior that we will discuss later. 

C. Extraction of phase shift 

Finally we turn to a major subject of this investigation, the extraction of the phase shift 
from the scattering process. Comparison with the result, eq. (0|, in the small amplitude 
approximation serves as a crucial test for the quality of the numerical solution to the time- 
dependent equation of motion. 

Numerically solving the equation of motion ()5j) to determine the full momentum depen- 
dence of the phase shift faces various obstacles. First we have to incorporate the above 
mentioned displacement of the kink in the integral, eq. ( ITT)) . This is straightforwardly ac- 
complished by restricting the integration interval to the regime of the wave-packet. For very 
late times (tf > 100) this regime is clearly separated from the kink. Other obstacles are more 
cumbersome. Wave-packet components with small momenta take a long time to finalize the 
interaction with the kink. Hence we need to solve the equation of motion for a large interval 
on the time axis. In the dimensionless units defined above we consider t £ [0,200]. We 
also vary the upper limit to ensure stability of the results. However, components with large 
momenta will propagate a long distance in the same time interval. Hence we also need to 
consider a large interval in coordinate space. In order to reliably find the Fourier transform, 
eq. (TT6|) we require a dense grid in coordinate space for large momentum components. This 
increases the numerical cost additionally. To keep the numerical effort within a manageable 
range, it is therefore appropriate to split the computation in (at least) two parts. To extract 
the phase shift for small momenta, we consider a large time interval but a small interval in 
coordinate space. This leads to unreliable results at large momenta. For that regime we 
consider a small time interval but a large one in coordinate space; together with a dense 
grid. At intermediate momenta the two procedures yield identical results. Furthermore we 
have the freedom to tune the parameters of the wave-packet, k and to suit the consid- 
ered regime in momentum space. These issues are indicated in figure In that figure we 
have matched two numerical treatments at k — 5. We also display the numerical result for 
the imaginary part of the right hand side of eq. ( TT6)) . Its deviation from unity serves as 
a further test on the numerical accuracy. As expected this occurs for very small and very 
large momenta. Otherwise the agreement with the analytical result, eq. (TT)) is astonishing. 
Certainly, further segmentation of the momentum axis and optimization in each segment 
will yield even better agreement. We have obtained the result displayed in figure El for a 
small amplitude (do = 0.05) for which the small amplitude approximation is expected to 
be accurate. Figure El shows that we find identical phase shifts for different amplitudes a^. 
Though we find differences in the imaginary part of the phase shift, we believe them to 
be short-comings of the numerical procedure and conclude that the scattering data do not 
exhibit non-linear effects even at moderate amplitudes. 

D. Non linear effects in the kink background 

As indicated in section [TTJ the extraction of the phase shift relies on the small-amplitude 
assumption. If Oirj^ ) effects are not negligible we cannot expect the right-hand-side of 
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FIG. 6: (Color online) Extraction of the phase shift for various momentum regimes. Left panel: 
do = 0.05, right panel: ag = 0.15. Different numerical parameters serve to improve the agreement 
with the analytical result in two distinct regimes. Left regime: optimized for small and moderate 
momenta, right regime: optimized for large momenta. 

eq. (I17|) to have unit magnitude. Rather we expect it to be less than one, corresponding to 
particle production. Yet we did not observe this effect even for moderate amplitudes, since 
the wave-packet remains compact, cf. figure For moderate amplitudes we have identified 
a small displacement of the kink as the single major effect of the non-linear dynamics. It 
corresponds to an attraction shortly before and a repulsion shortly after the interaction with 
the wave-packet. 

This displacement also complicates the extraction of that part of the energy density, 
e(x) that is associated with the wave-packet. Once this is properly done, e(x) serves as a 
probability distribution of the wave-packet. The propagation of its center does not exhibit 
consequences of the non-linear dynamics either. This propagation is not significantly altered 
by the presence of the kink. However, the spread of the wave-function shows some increase 
with the amplitude of the wave-packet. 

We show the typical behavior of the energy density, e(x, t) in figure[7]as we further increase 
the amplitude. Surprisingly, there is no footprint from the kink at Xq. Instead we observe 
that the energy density splits into two pieces of different velocities and there is a dominant 
peak in e(x, t) at the back of the more quickly propagating piece. The field configuration 
itself reveals the answer to this peculiarity. First we observe that also the wave-packet splits 
in (at least) two pieces of different velocities. This is an indication for particle production. 
Actually we also observe such splitting for the trivial background, however, it sets in at 
somewhat larger amplitude 5 . More notably, the transition from one vacuum configuration 
(0o = — 1) to the other (0o — 1), which is required to occur somewhere by conservation of 
the topological charge, now emerges at the back of the dominant piece of the wave-packet 
rather than at the position of the kink before the interaction. Obviously the dominant 
peak in e(x,t) results from the kink being dragged by the wave-packet. As a pronounced 
non-linear effect we find the transition from a stationary kink to a co-moving kink as the 
amplitude exceeds a certain value a c . We extract this value from the behavior of the energy 



5 For the vacuum background we did not observe it for clq = 0.25 but for oq = 0.35. 
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FIG. 7: (Color online) Results for real initial conditions using oo = 0.25. The numbers next to 
the lines refer to the time variable. Note that the second structure at t = 150 is partially covered 
by the main structure at t = 100. Left panel: subtracted energy density e(x,t), right panel: field 
configuration (p(x,t). 



density. When a structure persists at a value of d « xq that varies only marginally in time 
the kink is considered to be stationary. Disappearance of this structure provides the critical 
value. For ko = 4 and = 2 we find a c = 0.201. Obviously, it is impossible to extract 
scattering data from these structures. This is even more the case as parts of the initial 
wave-packet now trail after the kink. 



IV. SUMMARY AND OUTLOOK 



We have performed numerical simulations in the kink-soliton model. In particular we 
have set up the initial conditions to investigate the interaction of the kink with a wave- 
packet. Within this approach we have reconstructed the phase shift known from potential 
scattering in the small amplitude approximation. By choosing a significantly wide initial 
spectral function (in momentum space) this turned possible by a single integration of the 
equation of motion in coordinate space. However, numerical accuracy can be improved 
by optimizing both, numerical parameters and initial conditions, to a special momentum 
regime. This technique makes no particular reference to the background configuration and 
can hence be employed to other cases for which no analytic result is available. Although our 
technique also captures non-linear effects we did not observe them for the phase shift in the 
regime of small and moderate amplitude of the wave-packet. In this regime the kink does 
not pick up any energy from the wave-packet. However, the kink gets slightly displaced 
as a result of the interaction. Surprisingly, the displacement is opposite to the propagation 
direction of the kink. This displacement occurs as attraction before and repulsion after the 
interaction. Upon further increase of the amplitude the kink indeed picks up kinetic energy 
and gets dragged by the wave-packet. This is our main result regarding the search for 
effects of the non-linear dynamics. We interpret this effect as a signal for the existence of 
a critical strength of the amplitude beyond which the attraction between the wave-packet 
and the kink inescapable. Once the kink co-moves with the wave-packet, the scattering 
process cannot be uniquely identified nor can a phase shift be extracted. Furthermore the 
wave-packet splits into distinct pieces for sufficiently large amplitudes. We interpret this 
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splitting as particle production. It is also a consequence of the non-linear dynamics because 
it is equally observed without the kink background. 

Eventually we want to use these numerical methods to further study the interaction of 
a kink-antikink system. The comparison with the above results should shed some light on 
whether and how crossing symmetry manifests itself in (topological) soliton models. This is a 
challenging problem because unlike in perturbative calculations it does not simply imply the 
rotation of Feynman diagrams but rather relates topologically distinct sectors that are not 
connected by smooth transformations. The numerical analysis will closely follow the st udy 



of refs. [13|, |24j. The kink-antikink system is known to exhibit distinct features [lOj, [25 
depending on whether or not the relative velocity exceeds the critical value v c = 0.2598. 
Below v c the kinetic energy is not large enough to overcome the kink-antikink attraction and 
a second interaction occurs. Above v c the kinetic energy overcomes the barrier that causes 
the attraction. It will be interesting and challenging to establish a connection between the 
critical value of the amplitude (a c ) for the kink to take up kinetic energy and the above 
mentioned critical velocity (v c ). The existence of such a connection seems possible because 
both scenarios reflect a critical strength for sufficient attraction to result from the non-linear 
interaction within the model. 

The kink induces two bound states, one with energy u\ =v / 3^ and another (half-bound 
state) with energy u\ = It will be interesting to see whether there are connections 

between these bound states and the co-moving kink. If so, this co-moving object would 
correspond to an excited state of the kink. 
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